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Although the population of luminous quasars rises and falls over a period of ~ 10 9 
I years, the typical lifetime of individual quasars is uncertain by several orders of magni- 

tude. We show that quasar clustering measurements can substantially narrow the range 
of possible lifetimes with the assumption that luminous quasars reside in the most mas- 



sive host halos. If quasars are long-lived, then they are rare phenomena that are highly 
biased with respect to the underlying dark matter, while if they are short-lived they 
reside in more typical halos that are less strongly clustered. For a given quasar lifetime, 



we calculate the minimum host halo mass by matching the observed space density of 
quasars, using the Press-Schechter approximation. We use the results of Mo & White 
to calculate the clustering of these halos, and hence of the quasars they contain, as a 
function of quasar lifetime. A lifetime of tg = 4 x 10 7 years, the e-folding timescale 
of an Eddington luminosity black hole with accretion efficiency e = 0.1, corresponds 
to a quasar correlation length ro ~ 10/i _1 Mpc in low-density cosmological models at 



z = 2 — 3; this value is consistent with current clustering measurements, but these have 
large uncertainties. High-precision clustering measurements from the 2dF and Sloan 
quasar surveys will test our key assumption of a tight correlation between quasar lumi- 
nosity and host halo mass, and if this assumption holds then they should determine £q 
to a factor of three or better. An accurate determination of the quasar lifetime will show 
whether supermassive black holes acquire most of their mass during high-luminosity ac- 
cretion, and it will show whether the black holes in the nuclei of typical nearby galaxies 
were once the central engines of high-luminosity quasars. 

Subject headings: galaxies: quasars - cosmology: dark matter, large-scale structure of 
the universe 
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1. Introduction 

Mounting evidence for the existence of supermassive black holes in the centers of nearby 
galaxies (recently reviewed by, e.g., Richstone et al. 1998) supports the long-standing hypothesis 
that quasars are powered by black hole accretion (e.g., Salpeter 1964; Zel'dovich & Novikov 1964; 
Lynden-Bell 1969). However, one of the most basic properties of quasars, the typical quasar lifetime 
tQ, remains uncertain by orders of magnitude. The physics of gravitational accretion and radiation 
pressure provides one natural timescale, the e-folding time t e = Mrh/M = 4 x 10 8 e Z years of a 
black hole accreting mass with a radiative efficiency e = L/Mc 2 and shining at a fraction I = L/Le 
of its Eddington luminosity (Salpeter 1964). But while e ~ 0.1 and I ~ 1 are plausible values for 
a quasar, it is possible that black holes accrete much of their mass while radiating at much lower 
efficiency, or at a small fraction of Le- The task of determining tQ must therefore be approached 
empirically. 

The observed evolution of the quasar luminosity function imposes a strong upper limit on tQ 
of about 10 9 years, since the whole quasar population rises and falls over roughly this interval (see, 
e.g., Osmer 1998). The lifetime of individual quasars could be much shorter than the lifetime of 
the quasar population, however, and lower limits of tQ ~ 10 5 years rest on indirect arguments, 
such as the requirement that quasars maintain their ionizing luminosity long enough to explain 
the proximity effect in the Lya forest (e.g., Bajtlik, Duncan, & Ostriker 1988; Bechtold 1994). A 
typical lifetime tQ ~ 10 9 years would imply that quasars are rare phenomena, arising in at most 
a small fraction of high-redshift galaxies. Conversely, a lifetime as low as tQ ~ 10 5 years would 
imply that quasars are quite common, suggesting that a large fraction of present-day galaxies went 
through a brief quasar phase in their youth. 

The comoving space density &(z) of active quasars at redshift z is proportional to tQiiH(z), 
where nn is the comoving space density of quasar hosts. "Demographic" studies of the local black 
hole population (e.g., Magorrian et al. 1998; Salucci et al. 1999; van der Marel 1999) have opened 
up one route to determining the typical quasar lifetime: counting the present-day descendants of 
the quasar central engines in order to estimate uh(z) and thus constrain tQ by matching &(z). 
Roughly speaking, the ubiquity of black holes in nearby galaxies suggests that quasars are common 
and that tg is likely in the range 10 6 - 10 (e.g., Richstone et al. 1998; Haehnelt, Natarajan, & Rees 
1998; Salucci et al. 1999). However, as Richstone et al. (1998) emphasize, the lifetime estimated 
in this way depends crucially on the way one links the mass of a present-day black hole to the 
luminosity of a high-redshift quasar, which in turn depends on assumptions about the growth of 
black hole masses since the quasar epoch via mergers or low-efficiency accretion. 

In this paper we propose an alternative route to the quasar lifetime, using measurements 
of high-redshift quasar clustering. The underlying idea goes back to the work of Kaiser (1984) 
and Bardeen et al. (1986): in models of structure formation based on gravitational instability 
of Gaussian primordial fluctuations, the rare, massive objects are highly biased tracers of the 
underlying mass distribution, while more common objects are less strongly biased. Therefore, a 
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longer quasar lifetime tQ should imply a more clustered quasar population, provided that luminous 
quasars reside in massive hosts. The specific calculations that we present in this paper use the 
Press-Schechter (1974; hereafter PS) approximation for the mass function of dark matter halos and 
the Mo & White (1996, hereafter MW) and Jing (1998) approximations for the bias of these halos 
as a function of mass. The path from clustering to quasar lifetime has its own uncertainties; in 
particular, our predictions for quasar clustering will rely on the assumption that the luminosity of 
a quasar during its active phase is a monotonically increasing function of the mass of its host dark 
matter halo. However, the assumptions in the clustering approach are at least very different from 
those in the black hole mass function approach, and they can be tested empirically by detailed 
studies of quasar clustering as a function of luminosity and redshift. 

Our theoretical model of quasar clustering follows a general trend in which the study of quasar 
activity is embedded in the broader context of galaxy formation and gravitational growth of struc- 
ture (e.g., Efstathiou & Rees 1988; Turner 1991; Haehnelt & Rees 1993; Katz et al. 1994; Haehnelt, 
Natarajan, & Rees 1998; Haiman & Loeb 1998; Monaco, Salucci, & Danese 2000; Kauffmann & 
Haehnelt 2000). This paper also continues a theme that is prominent in recent work on the clus- 
tering of Lyman-break galaxies, namely that the clustering of high-redshift objects is a good tool 
for understanding the physics of their formation and evolution (e.g., Adelberger et al. 1998; Katz, 
Hernquist, & Weinberg 1999; Kolatt et al. 1999; Mo, Mao, & White 1999). Our model of the quasar 
population is idealized, but by focusing on a simple calculation with clearly defined predictions, we 
hope to highlight the link between quasar lifetime and clustering strength. After presenting the 
theoretical results, we will draw some inferences from existing estimates of the quasar correlation 
length. However, our study is motivated mainly by the anticipation of vastly improved measure- 
ments of quasar clustering from the 2dF and Sloan quasar surveys (see, e.g., Boyle et al. 1999; 
Fan et al. 1999; York et al. 2000). These measurements can test various hypotheses about the 
origin of quasar activity, including our primary assumption of a monotonic relation between quasar 
luminosity and host halo mass. If this assumption proves valid, then the first major physical result 
to emerge from the 2dF and Sloan measurements of high-redshift quasar clustering will be a new 
determination of the typical quasar lifetime. 



2. Method 
2.1. Overview 

We adopt a simple model of the high-redshift quasar population that is, doubtless, idealized, 
but which should be reasonably accurate for our purpose of computing clustering strength as a 
function of quasar lifetime. We assume that all quasars reside in dark matter halos and that a 
given halo hosts at most one active quasar at a time. The first assumption is highly probable, since 
a dark matter collapse is necessary to seed the growth of a black hole, and the second should be a 
fair approximation at high redshift, where the masses of large halos are comparable to the masses 
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of individual galaxy halos today. 

Our strongest and most important assumption is that the luminosity of a quasar during its 
active phase is monotonically related to the mass of its host dark matter halo, and that all suf- 
ficiently massive halos host an active quasar at some point. More precisely, we assume that an 
absolute-magnitude limited sample of quasars at redshift z samples the most massive halos present 
at that redshift, and that the probability that a halo above the minimum host mass M m \ n harbors 
an active quasar at any given time is tQ/tn, where tQ is the average quasar lifetime and tn is the 
halo lifetime. We can therefore compute the value of M m i n for a quasar population with comoving 
space density $(z) from the condition 



We compute n(M) using the PS approximation, and we compute the bias of halos with M > M m - m 
using the MW approximation. 

A connection between quasar luminosity and host halo mass is plausible on theoretical grounds 
— the cores of massive halos collapse early, giving black holes time to grow, and these halos provide 
larger gas supplies for fueling activity. It is also plausible on empirical grounds — local black hole 
masses are correlated with the host spheroid luminosity (Magorrian et al. 1998; van der Marel 
1999; Salucci et al. 1999), which in turn is correlated with stellar velocity dispersion (Faber & 
Jackson 1976). A precisely monotonic relation is certainly an idealization, and we explore the 
effects of relaxing this assumption in §4.1. The assumption of an approximately monotonic relation 
can be tested empirically by searching for the predicted relation between clustering strength and 
luminosity, as we discuss in §4.3. 

The ubiquity of black holes in luminous local spheroids supports our assumption that all 
sufficiently massive halos go through a quasar phase. However, once the quasar space density 
declines at z < 2, the occurrence of quasar activity must be determined by fueling rather than 
by the mere existence of a massive black hole, so it is not plausible that all large halos host a 
low-redshift quasar. We therefore apply our model only to the high-redshift quasar population, at 



We implicitly assume that a quasar turns on at a random point in the life of its host halo. In 
this sense, our model differs subtly from that of Haehnelt, Natarajan, & Rees (1998), who assume 
that a quasar turns on when the halo is formed, but this difference is unlikely to have a significant 
effect on the predicted clustering. Haehnelt, Natarajan, & Rees (1998) pointed out that a longer 
quasar lifetime would correspond to stronger quasar clustering because of the association with rarer 
peaks of the mass distribution, but they did not calculate this relation in detail. 

Because the quasar lifetime enters our calculation only through the probability tq/tn that 
a halo hosts an active quasar at a given time, it makes no difference whether the quasar shines 
continuously or turns on and off repeatedly with a short duty cycle (as argued recently by Ciotti 
& Ostriker 1999). For our purposes, tQ is the total time that the quasar shines at close to its peak 




(1) 



z > 2. 
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luminosity. We also assume that quasars radiate isotropically, with a beaming factor /# = 1, but 
because a smaller beaming factor simply changes the conversion between observed surface density 
and intrinsic comoving space density, all of our results can be scaled to smaller average beaming 
factors by replacing ig with fstQ. 



2.2. Notation 

All of our calculations assume Gaussian primordial fluctuations. We denote by P(k) the power 
spectrum of these fluctuations as extrapolated to the present day (z = 0) by linear theory. The 
rms fluctuation of the linear density field on mass scale M is 



a(M) = -— / dkk 2 P{k)W 2 {kr) 
2tt z Jo 

where 



1/2 



(2) 



— . 3( kr sin kr — cos kr) ( 3M \ ^ 3 . , 

W{kr) = (kr? ' r= U^J (3) 

is the Fourier transform of a spherical top hat containing average mass M. The mean density of 
the universe at z = is p = 2.78 x 10 n Q M h 2 M Mpc~ 3 , with h = H /(100 km s" 1 Mpc -1 ). The 
rms fluctuation can be considered as a function of either the mass scale M or the equivalent radius 
r. We define the normalization of the power spectrum by the value of a$ = a(r = 8/i _1 Mpc). 

The rms fluctuation of the linear density field at redshift z is 

a(M,z) = <t(M)D(z), (4) 

where D(z) is the linear growth factor D(z), defined so that D(z = 0) = 1. The general expression 
for the growth factor in terms of the scaled expansion factor y = (1 + z)^ 1 is 

where D(y) = 5{y) for Om = 1, D(y) = S(y)/S(l) for Qm < 1, and the dimensionless time variable 
is r = Hot (Heath 1977; Carroll, Press, & Turner 1992). If the dominant energy components are 
pressureless matter and a cosmological constant with ft a = A/3i?Q, then the Friedmann equation 
implies 



(^j =l + M (2 / - 1 -l)+n A (y 2 -l). 



(6) 



For an Om = 1, Qa = universe, D(z) = (1 + z)~ x . Peebles (1980, eq. 11.16) gives an exact 
analytic expression for D(z) for the case < 1, = 0, and Carroll, Press, & Turner (1992, 
eq. 29) give an accurate analytic approximation for ^ 0. In our notation, a(M) without any 
explicit z always refers to the rms linear mass fluctuation on scale M at z = 0. 
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At any redshift, we can define a characteristic mass M*(z) by the condition 

*[M m {z)]=6 e {z) = ^L, (7) 

where S c (z) is the threshold density for collapse of a homogeneous spherical perturbation at redshift 
z. Because we implicitly define the density field as "existing" at z = 0, the collapse threshold 5 c (z) 
increases with increasing redshift. For an £Im = 1, = universe, 

<y c ,o = 0.15(12tt) 2 / 3 » 1.69 

(see, e.g., Peebles 1980, §19). For other models, we incorporate the dependence of 5 Ct o on VLm in 
Appendix A of Navarro, Frenk, & White (1997), but because fij^f approaches one at high redshift 
in all models this correction to 5 C is less than 2% in all of the cases that we consider. 



2.3. From the Quasar Lifetime to the Minimum Halo Mass 



For a specified quasar lifetime, we compute the minimum halo mass by matching the comoving 
number density 3>(z) of observed quasars, accounting for the fact that only a fraction tQ/tn of host 
halos will have an active quasar at the time of observation. The matching condition is equation (1), 
or, putting in explicit mass and redshift dependences, 



[•CC 

Jm„ 



dM 



tQ 



t H (M,z) 



n(M, z) 



(8) 



If tQ > tji(M,z), we set the factor tQ/tn to unity. For the halo number density we use the PS 
approximation, 



n(M, z)dM = - 



2 P0 5c(z) da{M ) 
7r M a 2 (M) dM 



exp 



S 2 J 



2a 2 (M) 



dM, 



(9) 



where po is the mean density of the universe at z = 0, a(M) is the rms fluctuation given by 
equation (2), and 5 c (z) is the critical density for collapse by redshift z. 

In a gravitational clustering model of structure formation, halos are constantly growing by 
accretion and mergers, so the definition of a "halo lifetime" is somewhat ambiguous. For £Im(z) ~ lj 
a typical halo survives for roughly a Hubble time before being incorporated into a substantially 
larger halo, since the age of the universe at redshift z is also the characteristic dynamical time 
of objects forming at that redshift. Thus, to a first approximation, one could simply substitute 
tn(M,z) = tjj(z) in equation (9). We can do somewhat better by using the extended Press- 
Schechter formalism (e.g., Bond et al. 1991; Lacey & Cole 1993) to calculate the average halo 
lifetime, thereby accounting for the dependence of tjj on the power spectrum shape and the halo 
mass. Structure grows more rapidly in a cosmology with a redder power spectrum, and more 
massive halos accrete mass more rapidly. 



Equation (2.22) of Lacey & Cole (1993) gives the probability that a halo of mass M\ existing 
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at time t\ will have been incorporated into a new halo of mass greater than M 2 by time t 2 : 



P(S < S 2 ,U2\Si,wi) 



1 (uji - 2oj 2 ) 



-exp 



2uj 2 {uji - uj 2 ) 



1 - erf 



S 2 (wi - 2uj 2 ) + Siuj 2 
^j2S\S 2 {Si — S 2 ) 



+ - 



1 - erf 



S±uj 2 — S 2 uj\ 
y / 2SiS 2 (Si - 



(10) 



where S± = a 2 (Mi), S 2 = a 2 (M 2 ), lo\ = S c (ti), and uj 2 = 6 c (t 2 ). In this equation, lo plays the 
role of the "time" variable, with uj 2 < oj\ corresponding to t 2 > t±, and S plays the role of the 
"mass" variable, with S 2 < S\ corresponding to M 2 > Mi. For a halo of mass M existing at time 
tu(z), we define the halo lifetime to be the median interval before such a halo is incorporated into 
a halo of mass 2M. Thus, tn(M,z) = is — tjj(z), where is is the time at which the probability 
in equation (10) equals 0.5, for S\ = a 2 (M) and S 2 = a 2 (2M). Clearly other plausible definitions 
of tH(M,z) are possible, and they would give answers different by factors of order unity. With 
our definition, a black hole that lights up repeatedly is considered the "same" quasar as long as 
the mass of its host halo remains the same within a factor of two. If the host merges into a much 
larger halo and the black hole lights up again, it is considered a "new" quasar. We show the halo 
lifetimes for different masses and power spectra when we discuss specific models below. 

For comoving space densities &(z), we adopt values based on the work of Boyle et al. (1990), 
Hewett, Foltz, & Chaffee (1993), and Warren, Hewett, & Osmer (1994). Since observations con- 
strain the number of objects per unit redshift per unit solid angle, the conversion to comoving space 
density depends on the values of the cosmological parameters. We provide the formulas for these 
conversions in the Appendix, and in Table Al we list our adopted values of <E> (z) and the surface 
densities of objects to which these space densities correspond. In general, &(z) represents the space 
density of quasars above some absolute magnitude, corresponding to a surface density above some 
apparent magnitude. In §4.3 we discuss how to scale our results to predict the clustering of samples 
with different measured surface densities. 



2.4. From Minimum Halo Mass to Clustering Length 



Halos with M > M* are clustered more strongly than the underlying distribution of mass. 
MW give an approximate formula, 



b(M,z) = l + ^- 

<J C ,0 



a 2 {M) 



(11) 



for the bias factor of halos of mass M at redshift z. On large scales, the ratio of rms fluctuations in 
halo number density to rms fluctuations in mass should be b(M, z). This formula is derived from an 
extended Press-Schechter analysis, and it agrees fairly well with the results of N-body simulations on 
scales where the rms mass fluctuations are less than unity. The MW formula becomes less accurate 
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for halos with M < M*, i.e., a(M) < 6 c (z), and Jing (1998) provides an empirical correction that 
fits the N-body results, 



b{M, z) = [l + 



«5 C)0 [a 2 (M) 



^feH ■ (12) 



where n c g = 3 — 6 (dlncr/eZlnM) is the effective index of the power spectrum on mass scale M. 

According to our model, the quasars at redshift z only reside in halos of mass M > M m j n . The 
effective bias of these host halos is the bias factor (12) weighted by the number density and lifetime 
of the corresponding halos: 

Because the halo number density drops steeply with increasing mass, the effective bias is usually 
only slightly larger than the bias factor at the minimum halo mass, 6(M m i n , z). 

As our measure of clustering amplitude, we use the radius r\ of a top hat sphere in which 
the rms fluctuation gq of quasar number counts (in excess of Poisson fluctuations) is unity. This 
quantity is similar to the correlation length ro at which the quasar correlation function £(r) is 
unity, but it can be more robustly constrained observationally because it does not require fitting 
the scale-dependence of £(r). For a power law correlation function £(r) = (r/ro) -18 , r\ ps IAvq. 
With our adopted approximation for the bias, r\ is determined implicitly by the condition 

<7Q(n, z) = 6 eff (M min , z)a{r 1 )D{z) = 1, (14) 

where &(ri) is the rms linear mass fluctuation at z = in spheres of radius r\. For a specified 
cosmology, mass power spectrum P(k), quasar lifetime £q, and comoving space density &(z), we 
determine r\ from equation (14), computing a(r) from equation (2), D(z) from equation (5), M m \ n 
from equation (8), and 6 e ff(M m i n ) from equations (12) and (13). 



2.5. Results for Power Law Power Spectra 



Models with a power law power spectrum, P(k) oc k n , provide a useful illustration of our 
methods, since many steps of the calculation can be done analytically. For such models, the 
dependence of rms fluctuation on mass is also a power law, 



a{M) 



a(M, z) 



<5 c ,o 



D(z) D(z) 



M 



-(3+n)/6 



= S c (z) 



M 



-(3+n)/6 



(15) 



where M*(z) is the characteristic non- linear mass defined by equation (7). With this substitution, 
the PS mass function can be expressed as a function of M* (z) and the dimensionless mass variable 
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x = M/M*{z). Integrating to obtain the comoving number density of objects with mass M > M, 
yields 

+ A (P0_ 
vr V 6 H M > 



N(M > M min ) 



ax x e exp 



1 n+3 

"2 X 3 



(16) 



For power law models with £Im = 1> the ratio of the halo lifetime tu(M,z) to the age of 
the universe at redshift z depends only on M/M*(z) and has no separate dependence on redshift. 
Figure 1 shows tn and the ratio tn/tu as a function of M/M* for power law models with n = 0, 
— 1, and —2. More massive halos tend to accrete mass more quickly and therefore have shorter 
median lifetimes. At a given value of M/M*, the halo lifetime is shorter for lower n because a 
greater amount of large scale power causes the typical mass scale of non-linear structure to grow 
more rapidly. Although the calculation of the median halo lifetime via equation (10) is moderately 
complicated, the median lifetime for large masses asymptotically approaches a constant value 



t H (M,z) 



2 (3+n)/2 _ 1 



tu(z), M > M m (z), M = 1, 



(17) 



(Lacey & Cole 1993). We will find below that the predicted masses of quasar host halos are indeed 
in this asymptotic regime for most plausible parameter choices. The halo lifetime is longer for 
£Im < 1 than for Qm = 1 because fluctuations grow more slowly in a low-density universe, but in- 
still asymptotically approaches a constant value. The dotted curves in Figure 1 illustrate the case 
n = —1, VLm = 0.3, S7a = 0. The cosmological parameters for all of our models with power law 
power spectra are summarized in Table 1. 



The power law scaling of the rms fluctuation amplitude, equation (15), allows the bias for- 
mula (12) to be written 



b(M, z) 



1 + 



M 



5 C)0 [M*(z) 



(3+n)/3 



- 1 




M 



2 [M*(z)_ 



-2(3+n)/3 \ (0-06-0.02n) 

+ 1 . (18) 



Note that the second factor is very close to one for M > M*. Figure 2 shows the bias and the 
corresponding number- weighted effective bias (eq. [13]) as a function of M m i n /M*. For M m i n > M*, 
the effective bias is only slightly larger than 6(M mm ), since the number density of halos declines 
rapidly with increasing M. As equation (18) shows, the bias depends more strongly on M for larger 
n. However, the exponentially falling tail of the mass function at high M/M* is much steeper for 
higher n, as one can see from equation (16). As a result, the bias at fixed comoving number density 
is higher for smaller n in the high M/M* regime (see Fig. 3 below). 

Under the (good) approximation that the halo lifetime is given by the asymptotic formula (17) 
in the mass range of interest, the halo lifetime can be moved outside of the integral (8) for the 
number density of active quasars. The implied quasar lifetime as a function of minimum halo mass 
is then 

t (M \- 

Q[ mhl) ~ N(M > M min y 



(19) 
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Fig. 1. — Halo lifetimes vs. M/M* for the power law cosmologies listed in Table 1. Upper panels show the halo 
lifetimes in Gyr for each model at z — 2, 3, and 4. Lower panels show the ratio of the halo lifetime to the age of the 
universe. This ratio is independent of redshift for the ftju = 1 models, but not for the h = -1,Qm = 0.3 model. 



where N(M > M m j n ) is given by equation (16). For the Qm = 0.3, n = —1 model, we also 
use the asymptotic value of tjj, though this is no longer given by equation (17). We use a P(k) 
normalization ag = 0.5 for the three £Im = 1 models and ag = 1.0 for the 0,m = 0.3 model, in 
approximate agreement with the constraint on ag and implied by the observed mass function 
of rich galaxy clusters (White, Efstathiou, & Frenk 1993; Eke, Cole, k Frenk 1996). 

Equation (19) implicitly determines M m \ n /M*{z) given tQ. The top panels of Figure 3 show 
M m i n /M*(z) as a function of tcj/tu for z = 2, 3 and 4 and a constant comoving space density 
&(z) = 10~ 6 h 3 Mpc™ 3 . For the SIm = 1 cases, where tQ/tu depends only on n and M m i n /M*, 
the redshift dependence of M m i n /M„, arises solely from the presence of M* in the number density 
formula (16). As M* increases with, decreasing redshift, the value of x m j n = M m i n /M* must decrease 
to keep $(z) constant. Smaller values of n lead to higher values of M^^/M* because of the gentler 
fall off of the mass function at large M/M*. The n = —2 curves become flat at the largest values 
of tq/tu because tQ begins to exceed the halo lifetime tn, implying that all halos above M m - m are 
occupied by quasars. The difference between the open and Qm = 1 curves for n = —1 reflects 
mainly the larger values of ag and D{z) in the open model, which lead to a lower value of po/M* 
in the mass function (16) and therefore require a lower value of M m i n /M* to compensate. 
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Fig. 2. — Bias vs. M/M* for the power law models with n — (solid), —1 (short-dashed), and —2 (long-dashed). 
Lower curves show 6(M m i n ) computed from equation (12), and upper curves show the number-weighted effective bias 
(eq. [13]). 

The middle panels of Figure 3 show the effective bias fr e ff(M m i n , z) for the power law models. 
As already remarked, the bias at fixed space density and tQ/tjj is higher for redder power spectra 
(smaller n) because of the much higher values of Mmin/M*, despite the partially counterbalancing 
effect of the stronger dependence of bias on mass at larger n. Physically, the higher bias for redder 
spectra reflects the greater influence of the large scale environment on the amplitude of small scale 
fluctuations. For a given model, the bias increases with increasing redshift, reflecting the increase 
in Mmin/M*; the change, however, is quite modest. 

The rms number count fluctuation on comoving scale r is 

/ r \-(3+n)/2 

°"o( r i z ) = b c s(M miQ , z)a(r, z) = 6 cff (M min , z)a 8 D(z) ( g^-i Mpc ) ■ ( 20 ) 

The quasar clustering length r\ is the scale on which this rms fluctuation amplitude is unity, 

n = S^Mpc x [& eff (M min , z)a s D(z)] 2/iz+n) . (21) 

The bottom panels of Figure 3 present the main result of this Section, the dependence of n on 
quasar lifetime for our four power law models at z = 2, 3 and 4. As anticipated, the quasar clustering 
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Fig. 3. — Halo mass, bias, and clustering length for the power law models, as a function of quasar lifetime. Top 
panels show the minimum mass M m i n required to obtain a space density $(z) = W~ e 'h 3 Mpc~ 3 for a given value of 
tci/tu at z — 2, 3 and 4. Middle panels show the corresponding effective bias 6 e ff for halos M m i n . Bottom panels show 
the clustering lengths n as a function of tct/tu- The clustering length is the radius of a top hat sphere in which rms 
number count fluctuations (in excess of Poisson) are unity. For a power law correlation function f (r) = (r/ro) -1 ' 8 , 
n ?s 1.4 ro. 

length shows a strong dependence on quasar lifetime. The relation between r\ and ig depends on 
the power spectrum index n, so the shape of the power spectrum must be known fairly well to 
determine ig from measurements of r±. The clustering at fixed tg/it/ is substantially stronger in 
the open n = — 1 model than in the Q,m = 1 model because the underlying mass distribution is 
more strongly clustered (larger as and D(z)). 

For a specified value of Hm, the cluster mass function imposes a reasonably tight constraint on 
the normalization as- It is nonetheless interesting to explore the sensitivity of the predicted quasar 
clustering to this normalization. More intuitive than the erg-dependence is the equivalent relation 
between the quasar clustering length and the corresponding clustering length of the underlying 
mass distribution at the same redshift, 

r lm = Sh^Mpc x [a 8 D(zf /{3+n) . (22) 

Figure 4 plots this relation at z = 3 for the four power law cosmologies and ig = t\j (top curve), 
O.ltu, O.Olif/, and O.OOliry (bottom curve), for values of ag ranging from 0.2 to 2.0. The values 
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Fig. 4. — The dependence of the quasar clustering length riQ on the mass clustering length ri m , for the four power 
law models at z — 3. In each case, the four lines show, from top to bottom, the lifetimes <q = tu, 0-ltu, O.Oltu, and 
O.OOltu- Results are computed for normalizations running from = 0.2 to as = 2.0. Open circles show n m and 
riQ for our standard choices of as, listed in Table 1. If bias were independent of ri m , the lines would parallel the 
diagonal of the box, which has a slope of 1.0. 

of ri m that correspond to the as values in Table 1 are marked with open circles. If the bias did 
not change with ri m , then the quasar clustering length r\ would grow in proportion to ri m , and 
the curves in Figure 4 would parallel the diagonal of the box, which has a slope of 1.0. However, 
increasing r\ m increases M* and therefore requires a lower value of M^^/M^ to match the quasar 
space density. The correspondingly lower bias partially compensates for the larger ri m , making the 
curves in Figure 4 shallower than the box diagonal. For n = and large tq/tu (the highest solid 
curve), the minimum mass M m j n lies far out on the tail of a steeply falling mass function. In this 
regime, a change in M* requires only a small change in M m i n /M* to compensate, so there is little 
change in b e g with n m , and the curves approach the n oc ri m lines that would apply for constant 
bias. A similar argument explains the steepening of all curves towards low ri m , where the small 
values of M* put the value of M Iimi further out on the tail of the mass function. 
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3. Results for Cold Dark Matter (CDM) Cosmologies 

The results of §2.5 confirm our initial contention that quasar clustering can provide a good 
diagnostic of the typical quasar lifetime. However, they show that the predicted clustering length 
also depends on the shape of the mass power spectrum and on the value of Qm, which influences 
the cluster normalization of as at z = and (together with determines the growth factor 
D{z). Accurate determination of tg from measurements of quasar clustering therefore requires 
reasonably good knowledge of the underlying cosmology. Fortunately, many lines of evidence now 
point towards a flat, low-density model based on inflation and cold dark matter (see, e.g., the 
review by Bahcall et al. 1999). In particular, recent studies of the power spectrum of the Lya 
forest imply that the matter power spectrum has the shape and amplitude predicted by COBE- 
and cluster-normalized CDM models with £Im ~ 0.4 at the redshifts and length scales relevant to 
the prediction of quasar clustering (Croft et al. 1999; Weinberg et al. 1999; McDonald et al. 2000; 
Phillips et al. 2000). 

For the power spectrum of our CDM models, we adopt P(k) oc k np T 2 (k) with scale-invariant 
{n p = 1) primeval inflationary fluctuations and the transfer function parameterization of Bardeen 
et al. (1986), 

T(k) = ln(1 9 + o!" 34g) t 1 + 3 - 89 ? + ( 16 - 1( ?) 2 + ( 5 - 46 9) 3 + (G.7lq) 4 } " V4 . (23) 

Here q = k/T and T, with units of (/i _1 Mpc) _1 , is the CDM shape parameter, given approximately 
by T = flMhexp(-Qb - \Z2hQb/QM) (Sugiyama 1995). We calculate a(M) and (da/dM) by 
numerical integration of this power spectrum. 

We consider five different CDM models with the parameters listed in Table 2. These models 
are chosen to illustrate a range of cosmological inputs and also to isolate the effects of different 
parameters on quasar clustering predictions. The rCDM, OCDM, and ACDM models have T = 0.2, 
in approximate agreement with the shape parameter estimated from galaxy surveys (e.g., Baugh 
& Efstathiou 1993; Peacock Sz Dodds 1994), and they have as values consistent with the cluster 
mass function constraints of Eke, Cole, & Frenk (1996). The tCDM and ACDM models arc 
approximately COBE-normalized. COBE normalization would imply a lower as for OCDM, but 
a slight increase in n p could raise ag without having a large impact on the shape of P(k) at the 
relevant scales. The OCDM and ACDM models are consistent with the Lya forest power spectrum 
measurements of Croft et al. (1999), but the rCDM model is not. OCDM is inconsistent with 
the observed location of the first acoustic peak in the cosmic microwave background anisotropy 
spectrum (e.g., Miller et al. 1999; Melchiorri et al. 1999; Tegmark &: Zaldarriaga 2000), and of the 
three models, only ACDM is consistent with the Hubble diagram of Type la supernovae (Riess et 
al. 1998; Perlmutter et al. 1999). 

We will use the comparison between the rCDM and SCDM models, with T = 0.2 and T = 0.5, 
respectively, to illustrate the impact of power spectrum shape for fixed Qm and as- The SCDM 
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Table 1. Power Law Model Parameters 
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Note. — Parameters of the four power law cosmological models discussed in §2.5. 
Columns 1-4 list the power spectrum index n, normalization <t 8 , scaled Hubble 
constant h, and mass density parameters Qm- Columns 5-7 list the values of 
(eq. [7]) at z = 2, 3, and 4, in units of h~ 1 M Q . 



Table 2. CDM Model Parameters 
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Note. - Parameters of the five CDM models discussed in §3. Column 1 lists 
the model name, columns 2-5 the power spectrum normalization and cosmologi- 
cal parameters, and column 6 the power spectrum shape parameter (see eq. [23]). 
Columns 7-9 list the values of M* (eq. [7]) at z = 2, 3, and 4, in units of h~ 1 M & . 
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model is cluster-normalized, but its as = 0.5 is well below the value as ~ 1.2 implied by COBE for 
n p = 1, r = 0.5 (e.g., Bunn & White 1997). The OCDM and ACDM models have the same P(k) 
shape and the same P(k) amplitude at z = 0, but at high redshift the OCDM model has stronger 
fluctuations because of a larger D(z). We therefore include the model ACDM2, which has as chosen 
to yield the same power spectrum amplitude as OCDM at z = 3. Differences between OCDM and 
ACDM2 isolate the impact of a cosmological constant for fixed high-redshift mass clustering. 

Figure 5 shows tn in Gyr (upper panels) and tii/tu (lower panels) as a function of M/M* for 
the CDM models at z = 2, 3, and 4. In contrast to the power law models shown in Figure 1, the 
ratio tff/tu does not approach a constant value but instead increases at very large M/M*. This 
increase can be understood with reference to the power law case: the effective power law index 
n e fr = 3 — 6 (d In a /din M) increases with increasing mass in a CDM spectrum, and larger values 
of n e fr correspond to slower growth of mass scales (and larger tn/tu) as shown in Figure 1. The 
difference between the SCDM and rCDM curves in Figure 5 reflects the higher values of n e g for 
the r = 0.5 power spectrum. The differences between the various T = 0.2 models largely reflect 
the differences in M*, and hence the differences in n e g at fixed M/M*, and they also reflect the 
differences in fluctuation growth rates. 
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Fig. 6. — Effective bias as a function of minimum halo mass, as in Figure 2, for the CDM models at z = 3. 

Figure 6 plots the effective bias against M m i n /M* for the five CDM models at z = 3. Figure 2 
showed that the value of b e g at fixed M m i n /M* is higher for larger n. The lines in Figure 6 curve 
upwards because n e g increases with mass scale, and to a good approximation the value of 6 c ff i n 
the CDM models equals the value of b e g at the same M m - m /M* in a power-law model of index 
n e ff(M m i n ). The difference between SCDM and rCDM in Figure 6 therefore reflects the higher n e fr 
values in SCDM, and the differences among the other models reflect the different values of M*, and 
hence the different values of n e s at fixed M m - m /M*. 

The top three panels of Figure 7 show the dependence of M m \ n /M* on ig at z = 2, 3, and 4; the 
values of are listed in Table 2. The calculation of M m j n via equation (8) incorporates both the 
dependence of halo lifetime on mass and the influence of and Q\ on the value of &(z) inferred 
from the quasar surface density (as discussed in the Appendix). The two CIm = 1 models have the 
lowest values of M* because of their lower ag and D(z), so they require the largest M m i n /M* to 
match the observed ®(z). The value of M* is smaller for ACDM than for OCDM because D(z) is 
smaller for the flat model, so ACDM requires larger M m i n /M*. The higher normalization of the 
ACDM2 model largely removes this difference, since a%D(z = 3) is matched to that of the OCDM 
model, but ACDM2 still has a slightly lower M* because of the influence of on S c (z). As a 
result, the M min /M* curve for ACDM2 lies just above that of OCDM at z = 3. 
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Fig. 7. — Minimum halo mass, effective bias, and clustering length as a function of £q, for the CDM models at 
z = 2, 3, and 4. Format is the same as Figure 3. 

The middle panels of Figure 7 show the effective bias values, which display the same relative 
dependence on tQ and cosmology as the M m i n /M* values. Because M m i n /M* > 1 in all of the CDM 
models, even for tQ as low as 10 5 years, the MW bias formula (11) yields nearly identical results 
to Jing's (1998) corrected formula (eq. [12]). 

The bottom panels of Figure 7 present the main results of this paper: the relation between the 
clustering length r\ and the quasar lifetime tQ for CDM models at z = 2, 3, and 4. The clustering 
length is an increasing function of quasar lifetime for the reasons outlined in the Introduction and 
detailed in §2. A longer tQ implies that quasar host halos are rarer, more highly biased objects. 
The change in the r\ vs. tQ relation with redshift reflects the evolution of the quasar space density 
and of the underlying mass fluctuations. For a given model and tQ, the predicted quasar clustering 
is weakest at z = 3, the peak of the quasar space density. The smaller quasar abundance at z = 4 
implies a higher bias of the host halo population, which more than compensates for the slightly 
weaker mass clustering. The clustering length grows between z = 3 and z = 2 because of both the 
drop in quasar space density and the growth of mass clustering. The T\ vs. tQ relation becomes 
flat at the largest tQ for the SCDM model at z = 4 and for the rCDM model at z = 3 and 4, where 
tQ exceeds the halo lifetime iff(M m i n ) and the value of M m \ D required to match $>(z) therefore 
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becomes independent of tQ. 

The differences between models reflect the differences in bias factors discussed above and the 
differences in the mass clustering. There are also differences in the values of inferred from the 
observed quasar surface density (see Appendix), but these have relatively little effect. The main 
separation in Figure 7 is between the low-density models and the fi^ = 1 models, which have 
weaker mass clustering because of their lower values of as and D(z). The VLm = 1 models have 
larger bias factors, but these are not enough to compensate for the smaller mass fluctuations. The 
r\ vs. tQ relations are also shallower for the Q,m = 1 models, because the values of M m j n lie further 
out on the steep, high-mass tail of the mass function, where a smaller change in M m - m can make 
up for the same change in tQ. The three low-density models yield very similar predictions. 

To facilitate comparison of future observational results to our predictions, we have fit polyno- 
mials of the form 



to each of the r\ vs. tQ curves shown in Figure 7. The values of the coefficients are given in 
Table 3, and the coefficients ao and a\ have the same value over the entire range in tQ. These fits 
are accurate to better than 3% in r\ for given tQ, or better than 10% in tg given n, for all cases 
except the SCDM and rCDM models at z = 4, where the maximum errors are 3% in r\ given tQ 
and 20% in tQ for given r\. 



As already mentioned in §2.3, the definition of a "halo lifetime" is somewhat ambiguous. We 
have so far adopted a definition of tu as the median time before a halo of mass M is incorporated 
into a halo of mass 2M. If we increase this mass ratio from 2 to 5 (a rather extreme value), then 
the typical halo lifetimes in our CDM models increase by factors of 2 — 4. Since it is the ratio tQ/tu 
that enters our determination of M m i n (eq. [8]), and hence fixes the bias factor, this change in tu 
would require an equal increase in tQ to maintain the same clustering length r\ . We conclude that 
the ambiguity in halo lifetime definition introduces a factor ~ 2 uncertainty in the determination 
of tQ from clustering measurements, in the context of our model. 

We have also assumed that quasar luminosity is perfectly correlated with host halo mass, so 
that matching the space density of an absolute-magnitude limited sample imposes a sharp cutoff 
in the host mass distribution at M = M m ; n . If there is some scatter in the luminosity-host mass 
relation, then some halos with M < M m j n will host a quasar above the absolute-magnitude limit 




(24) 



4. 



Discussion 



4.1. Sensitivity to model details 



-20- 



and some halos with M > M m i n will not. We can model such an effect by introducing a soft cutoff 
into equation (8): 

Hz) = f~ d M 3 (M) J -{SL_n(M,z) (25) 



with 

g(M) 



for M < ^ 



{ m^-i) ) M ~ ah for^<M<aM min (26) 
1 for M > aM m]n 



and a > 1. Adopting a soft cutoff slightly decreases M m \ n and, more significantly, reduces the 
value of 6 e ff by allowing some quasars to reside in lower mass halos, which are less strongly biased. 
Quantitatively, we find that setting a = 2, which corresponds to including halos down to M = 
M m - m /2, decreases the clustering length by < 6% for the shortest quasar lifetimes and < 10% for 
the longest quasar lifetimes. Matching a fixed n with an a = 2 cutoff requires lifetimes that are 
longer by a factor ~ 1 — 1.5 at short tQ and ~ 2 — 2.5 at long tQ. Longer lifetimes are more sensitive 
to scatter in the luminosity-host mass relation because b e g depends more strongly on M m i n /M* for 
these rarer objects. The assumption of a perfectly monotonic relation between quasar luminosity 
and host mass leads to the smallest tQ for a given r±. Thus if any scatter does exist in this relation, 
our model predictions for tQ effectively become lower limits to the quasar lifetime. 

Another simplification of our model is the assumption that a quasar is either "on" or "off" - 
each quasar shines at luminosity L for time tQ, perhaps divided among several episodes of activity, 
and the rest of the time it is too faint to appear in a luminous quasar sample. More realistically, 
variations in the accretion rate and radiative efficiency will cause the quasar luminosity to vary, 
especially if the black hole mass itself grows significantly during the luminous phase. Nonetheless, 
the maximum luminosity will still depend on the maximum black hole mass. At a given time, the 
luminous quasar population will include black holes shining at close to their maximum luminosity 
and "faded" black holes of higher mass. Because the host halos lie on the steeply falling tail of 
the mass function, the first component of the population always dominates over the second, and 
we therefore expect our clustering method to yield the time tQ for which a quasar shines within a 
factor ~ 2 of its peak luminosity. More strongly faded quasars are too rare to make much difference 
to the space density or effective bias. 

To illustrate this point, we consider the model of Haehnelt, Natarajan, & Rees (1998) in which 
a quasar hosted by a halo of mass M has a luminosity history L(t) = L${M) exp(— t/tQ), with 
a maximum luminosity L$(M) = aM proportional to the halo mass. In this model, the time 
that a quasar shines above the luminosity threshold L mm = Lo(M m i n ) of a survey is the visibility 
time t f Q = tQ ln(M/M m i n ). We can calculate M m - m for a given space density by substituting t'g for 
tQ in equation (8), then calculate 6 e fr(-Mmm) by multiplying the integrands in the numerator and 
denominator of equation (13) by the visibility weighting factor ln(M/M m j n ). The middle curves in 
Figure 8 compare r\{to) for the on-off {solid line) and exponential decay {dotted line) models, in 
the case of ACDM at z = 3 with our standard $(z). The curves are remarkably similar, showing 
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Fig. 8. — Clustering length vs. tQ for the ACDM model at z — 3 for two different models of quasar luminosity 
evolution at three different space densities &(z). The "on - off" model (solid lines) assumes the quasar luminosity is 
constant throughout its lifetime tQ and is the standard model we discuss in this paper. The central line shows results 
for the ACDM model at z — 3 with our standard <&(z). The other solid curves, related to the first by simple horizontal 
shifts, show results for space densities different by factors of 10 and 1/10 (bottom and top). In the exponential model 
(dotted lines), the quasar luminosity starts at some maximum luminosity proportional to the halo mass and decays 
with an e- folding timescale tQ. The middle line again corresponds to our standard $(2), and the other two dotted 
curves show results for space densities different by factors of 10 and 1/10 (bottom and top). 

that the lifetime inferred from clustering assuming an on-off model would be close to the e-folding 
timescale in an exponential decay model. The curves for the exponential decay model are slightly 
shallower because at low M m - m (low tQ) the mass function is not as steep, allowing faded quasars in 
more massive halos to make a larger contribution to b e g and thereby raise r±. Although results for 
a different functional form of L(M, t) would differ in detail, we would expect the lifetime inferred 
from clustering to be close to the "half-maximum" width of the typical luminosity history, for the 
general reasons discussed above. 

As mentioned in §2.1, we assume that quasars radiate isotropically. If they radiate instead 
with an average beaming factor /b < 1, then the true value of $(z) is larger than the observed 
value by a factor f^ 1 . The implied lifetime for a given n would therefore be larger by a factor f^ 1 
as well. 
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4.2. Interpretation of Existing Data 

After several attempts (Osmer 1981; Webster 1982), quasar clustering was first detected by 
Shaver (1984), and later by Shanks et al. (1987) and Iovino k, Shaver (1988). However, measure- 
ments of quasar clustering are still hampered by small, sparse samples, and even the best studies 
to date yield detections with only several-c significance. Given the limitations of current data, it is 
not surprising that different authors reach different conclusions about the strength of clustering and 
its evolution. Analyzing a combined sample of quasars with 0.3 < z < 2.2 from the Durham/AAT 
UVX Survey, the CFHT survey, and the Large Bright Quasar Sample, Shanks & Boyle (1994) 
and Croom & Shanks (1996) find a reasonable fit to the data with an VLm = 1 model that has 
£(r) = (r/ro)~ 7 , 7 ~ 1.8, and a constant comoving correlation length tq = 6/i _1 Mpc. La Franca 
et al. (1998) report a higher correlation length, r$ = 9.1 ± 2.0/i~ 1 Mpc for a 7 = 1.8 power law, in 
their 1.4 < z < 2.2 sample. 

If we adopt ro « 8/i _1 Mpc at z = 2 and a corresponding 77 w ll/i _1 Mpc, then the implied 
quasar lifetime is ~ 10 7 ' 5 years for the rCDM model and ~ 10 8 years for SCDM. The ro values 
quoted above are for Qm = 1> and because quasar pair separations are measured in angle and 
redshift, they should be increased by a factor ~ 1.5 in an Qm = 0.3, $7 a = 0.7 universe and a factor 
~ 1.3 in an Qm = 0.3, Oa = universe (roughly the inverse cube-roots of the volume ratios listed 
in Table Al). Adopting 77 ps 16/i _1 Mpc implies a lifetime tQ ~ 10 7 ~ 7 ' 5 years in our low-density 
models. However, these numbers must be considered highly uncertain because of the limitations of 
current data and because the space densities of the various observational samples do not necessarily 
match those assumed in our model predictions. 

All of these measurements are based mainly on quasars with z < 2. At higher redshift, Kundic 
(1997) and Stephens et al. (1997) have investigated clustering in the Palomar Transit Grism Survey 
(PTGS; Schneider, Schmidt, & Gunn 1994). Fitting a 7 = 1.8 power law, Stephens et al. (1997) 
find ro = 17.5 ± 7.5/i _1 Mpc for z > 2.7. This high correlation length (inferred from the presence of 
three close pairs in a sample of 90 quasars) could be a statistical fluke, but in the context of our 
model it is tempting to see it as a consequence of the high luminosity threshold of the PTGS survey, 
which might lead it to pick out the most strongly clustered members of the quasar population. 

4.3. Prospects 

The 2dF (Boyle et al. 2000; Shanks et al. 2000) and Sloan (York et al. 2000) quasar surveys 
will transform the study of quasar clustering over the next several years, yielding high-precision 
measurements for a wide range of redshifts. These measurements will allow good determination of 
the typical quasar lifetime tQ in the context of the model presented here. They will also test the key 
assumption of this model, the monotonic relation between quasar luminosity and host halo mass, 
by characterizing the clustering as a function of redshift and, especially, as a function of quasar 
absolute magnitude. 
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Figure 8 illustrates this test for the ACDM model at z = 3. Brighter quasars have a lower 
space density &(z), so they should have a higher minimum host halo mass M m [ n , and, because 
of the higher bias of more massive halos, they should exhibit stronger clustering. Fainter, more 
numerous quasars should exhibit weaker clustering. Figure 8 shows the predicted r\ vs. ig relation 
for samples with 1/10 and 10 times the space density of our standard case (3.42 quasars per square 
degree per unit redshift; see Table Al). In our standard on-off model {solid lines), a change in <E>(z) 
in equation (8) can be exactly compensated by changing ig by the same factor, so the solid curves 
in Figure 8 are simply shifted horizontally relative to each other. Our predictions in Figure 7 (see 
eq. [24]) can therefore be transformed to any quasar space density by changing tn in proportion to 
<3?(z). In the exponential decay model (dotted lines), the scaling of ig with &(z) is no longer exact, 
though it is still a good approximation. 

If there is a large dispersion in the relation between quasar luminosity and host halo mass, 
then the dependence of clustering strength on quasar space density will be much weaker than 
Figure 8 predicts. Detection of the predicted trend between luminosity and clustering, or definitive 
demonstration of its absence, would itself provide an important insight into the nature of quasar 
host halos. More generally, the parameters of a model that incorporates scatter (such as the a 
prescription of equation [26]) could be determined by matching the observed relation between n 
and ®(z). 

If the observations do support a tight correlation between luminosity and host halo mass, then 
the first property of quasars to emerge from the 2dF and Sloan clustering studies will be the typical 
lifetime tQ. For the low-density models in Figure 7, the slope of the correlation between r\ and 
log 10 tg is ~ 10, so a determination of r\ with a precision of 2/i _1 Mpc would constrain tQ to a 
factor of 10 - 2 « 1.6, for a specified cosmology. By the time these quasar surveys are complete, 
a variety of observations may have constrained cosmological parameters to the point that they 
contribute negligible uncertainty to this constraint. Instead, the uncertainty in tQ will probably be 
dominated by the limitations of the quasar population model, e.g., the approximate nature of the 
assumptions that the quasar luminosity tracks the halo mass, that there is only quasar per halo, and 
that the average lifetime tQ is independent of quasar luminosity. These assumptions can be tested 
empirically to some degree, but not perfectly. Despite these limitations, it seems realistic to hope 
that tQ can be constrained to a factor three or better by high-precision clustering measurements, 
a vast improvement over the current situation. It is worth reiterating that our assumption of a 
perfectly monotonic relation between luminosity and halo mass leads to the smallest tQ for an 
observed r±, since with a shorter lifetime there are simply not enough massive, highly biased halos 
to host the quasar population. 

A determination of ig to a factor of three will be sufficient to address fundamental issues about 
the physics of quasars and galactic nuclei. Comparison of ig to the Salpeter timescale will answer 
one of the most basic questions about supermassive black holes: do they shine as they grow? If 
>4x 10 7 years, the e-folding timescale for L ~ Le, e ~ 0.1, then quasar black holes increase 
their mass by a substantial factor during their optically bright phase. If ig is much shorter than 
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this, then the black holes must accrete most of their mass at low efficiency, or while shining at 
L Le- A short lifetime could indicate an important role for advection dominated accretion 
(Narayan, Mahadevan, & Quatert (1998) and references therein), or it could indicate that black 
holes acquire much of their mass through mergers with other black holes, emitting binding energy 
in the form of gravitational waves rather than electromagnetic waves. A determination of tq would 
also resolve the question of whether the black holes in the nuclei of local galaxies are the remnants 
of dead quasars. For example, Richstone et al. (1998) infer a lifetime tQ ~ 10 6 years by matching 
the space density of local spheroids that host black holes of mass M > 4 x 1O 8 M to the space 
density of high-redshift quasars of luminosity L^(M) > 6 x 10 46 erg s" 1 . If clustering implies a 
much longer lifetime, then these numerous local black holes may once have powered active nuclei, 
but they were not the engines of the luminous, rare quasars. 

We have assumed in our model that quasar activity is a random event in the life of the 
parent halo. Quasar activity might instead be triggered by a major merger, by a weaker "fly-by" 
interaction, or by the first burst of star formation in the host galaxy. Regardless of the trigger 
mechanism, the lifetime will be the dominant factor in determining the strength of high-redshift 
quasar clustering, if our assumed link between luminosity and halo mass holds. However, different 
triggering mechanisms might be diagnosed by more subtle clustering properties, such as features in 
the correlation function at small separations, or higher-order correlations. At low redshift, where 
the evolution of the quasar population is driven by fueling rather than by black hole growth, the 
nature of the triggering mechanism might play a major role in determining quasars' clustering 
properties. The calculations presented here illustrate the promise of quasar clustering as a tool for 
testing ideas about quasar physics, a promise that should be fulfilled by the large quasar surveys 
now underway. 

We thank James Bullock, Alberto Conti, Jordi Miralda-Escude, Patrick Osmer, and Simon 
White for helpful discussions. We also thank the referee for constructive suggestions, which led to 
our consideration of the exponential decay model in § 4. As we were nearing completion of this 
work, we learned of a similar, independent study by Z. Haiman & L. Hui (Haiman & Hui 2000); 
our general conclusions are consistent with theirs, although the approaches are quite different in 
detail, precluding a precise comparison of results. This work was supported in part by NSF grant 
AST96-16822 and NASA grant NAG5-3525. 



A. Converting from observed quasar numbers to $(z) 

The observed quantity that is measured in studies of quasar clustering and the quasar space 
density is the number of sources brighter than a given apparent magnitude m per unit redshift per 
unit solid angle on the sky. This surface density per unit redshift can be converted into a comoving 
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space density of objects brighter than a given absolute magnitude M, 

dN(< m) dQdz . A . 

^ <M ^^mrwM' (A1) 

where dV c (z) is the differential comoving volume element corresponding to dQ dz. Following the 
notation in Hogg (1999), this volume element is 



dV c (z) = D H ^f-dQdz, 
where Dh = c/Hq is the Hubble distance, Dm is the transverse comoving distance, 



(A2) 



D h sinhlv 7 ^^] for Q k > 
DhJ*^ iorQ k = 

D H ^^ S m[y^h\^] forO fe <0 



(A3) 



|«fc| 



and E(z) = y / ^m(1 + z) 3 + Q k (l + z) 2 + Q A , where Q k = 1 — Q M — ^A- For Q M = 1, = 0, the 
differential comoving volume element is 



dVc(z) = 4 
per steradian per unit redshift. 



c 



(1 + *) 



-3/2 



dQ dz 



(A4) 



The fact that dV c (z) depends on the cosmological parameters means that a given measured sur- 
face density of sources corresponds to a different comoving space density for different cosmological 
model parameters. The space density of quasars is commonly quoted for an Qm = 1 universe. To 
convert this space density (in units of h 3 Mpc~ 3 ) into the space density for a model with different 
values of Qm and Qa requires a correction of the form 

f(z, ^mi ^a) 



$(z,q m ,Qa) = <f>(z,n' M ,n' A ) 



/(2,Hm,0a) : 



where 



f(z, q m ,Qa) 



DhD 2 m 
E(z) ' 



(A5) 



(A6) 



This procedure converts the reported space density under one set of cosmological parameters back 
into the observed surface density and then converts the surface density into the space density for the 
new set of cosmological parameters. In the notation of Popowski et al. (1998), f(z, Qm, Qa) = gxf 2 , 
where / and g are given by their equations (5) and (6), respectively. 



In Table Al we list the factors to convert the space density in column 2, which is listed for 
Qm = 1, ^A = 0, to the corresponding space densities for Qm = 0.3, Qa = 0.0 and Qm = 0.3, Qa = 
0.7. The factors in Table Al are all less than unity because the comoving volume element is smallest 
in an Qm = 1 universe. 
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Table 3. r\ vs. tQ fitting coefficients 



Model 


a 


ai 


a 2 


a 


ai 


a>2 


a 


ai 


02 






z = 2 






z = 3 






z = 4 




SCDM 


12.47 


2.437 


0.1001 


10.62 


1.932 


0.0601 


11.26 


1.590 


-0.0042 


rCDM 


16.54 


3.163 


0.1662 


13.89 


2.422 


0.0882 


14.48 


1.964 


-0.0166 


OCDM 


27.18 


7.002 


0.7798 


25.32 


6.029 


0.5347 


30.26 


5.457 


0.1640 


ACDM 


26.79 


6.173 


0.4913 


23.48 


4.974 


0.3122 


26.12 


4.248 


0.0823 


ACDM2 


29.84 


7.441 


0.7465 


26.51 


6.093 


0.4840 


30.23 


5.274 


0.1393 



Note. Coefficients for polynomial fits (eq. [24]) to the predicted relations between 
quasar lifetime and clustering length shown in Figure 7. Column 1 lists the model name, 
columns 2-4 the coefficients for z = 2, columns 5-7 the coefficients for z = 3, and columns 
8-10 the coefficients for z = 4. 



Table Al. Quasar Space and Surface Density 



z 


*(z) 


dN 
dz dQ 


/(i.o.o.o) 

/(0.3.0.0) 


/(1.0.0.0) 
/(0.3.0.7) 


2 


1.889 • 10~ 6 


2.132 


0.425 


0.279 


3 


3.331 • 10~ 6 


3.417 


0.339 


0.253 


4 


3.200 • 10~ 7 


0.287 


0.287 


0.241 



Note. — Adopted values of the space density of quasars at z = 2, 3, and 
4, and cosmological conversion factors. Values of &(z) in column 2 are from 
Warren, Hewett, & Osmer (1994) (assuming £Im = 1) for quasars with M c < 
—24.5 (absolute continuum flux at 1216 A), converted from their adopted 
h = 0.75 to /i 3 Mpc~ 3 . Column 3 lists the abundance of quasars in number 
per square degree per unit redshift to which the space density in column 2 
corresponds. Columns 4 and 5 contain the ratios of the factors /(Ami^a) 
defined in equation (A6) needed to convert the space density in column 2, 
which is valid for Qm = 1, to space densities for the OCDM and ACDM 
models, respectively. 
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